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Abstract. The Brackett decrement in the Class I source HHIOO IR has been observed and analyzed to set con- 
straints on the origin of the IR HI emission in this young object. We have used both low resolution (R~800) 
observations of the Brackett lines from Br7 to Br24, and medium resolution (R~9000) spectra of the Br7, Brl2 
and Brl3 lines. The dereddened fluxes indicates that the lines remain moderately thick up to high quantum num- 
bers. Moreover, the profiles of the three lines observed in medium resolution are all broad and nearly symmetric, 
with a trend for the lines at high n-number to be narrower than the Br7 line. With the assumption that the 
three lines have different optical depths and consequently trace zones at different physical depths, we interprete 
the observed profiles as evidence that the ionized gas velocity in the HI emitting region is increasing as we move 
outwards, as expected in an accelerating wind more than in an infalling gas. 

We have modelled the observed line ratios and velocities with a simplified model for the HI excitation from a 
circumstellar gas with a velocity law V = Vq + {Vmax — Vq){1 — (ri/r)"). Such a comparison indicates that the 
observations are consistent with the emission coming from a very compact region of 4-6 R©, where the gas has 
been already accelerated to velocities of the order of 200 kms~^ , with an associated mass fiow rate of the ionized 
component of the order of 10~^ M0yr~^. This implies that the observed lines should originate either from a 
stellar wind or from the inner part of a disk wind, providing that the disk inner truncation radius is close to the 
stellar surface. It is also expected that the gas ionization fraction is relatively high as testified by the high rate of 
ionized mass loss derived. Our analysis, however, does not resolve the problem of how to reproduce the observed 
symmetrical line profiles, which at present are apparently difficult to model by both wind and accretion models. 
This probably points to the fact that the real situation is more complicated than described in the simple model 
presented here. 
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1. Introduction 

The emission from hydrogen recombination lines repre- 
sents the most direct manifestation of circumstellar ac- 
tivity in young stars. Lines from the Balmer series have 
been the main spectroscopic tool for identifying classi- 
cal T Tauri and Herbig AeBe stars. In spite of the large 
amount of observational data, however, the HI lines are 
so easily excited that it is very difficult to clearly define 
their origin. Originally, they were interpreted as being ex- 
cited in circumstellar ionized winds, on the basis of the 
often observed P Cygni profile exhibited by the Ha line 
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(Hartmann et al. 1990, Calvet et al. 1992). More recently, 
such an interpretation has been challanged due to the va- 
riety of the observed Balmer line profiles (e.g. Edwards et 
al. 1994, Reipurth et al. 1996) that are difficult to repro- 
duce by means of wind models, and the HI emission has 
been interpreted also in the framework of magnetospheric 
accretion models (Calvet & Hartman 1992, MuzeroUe et 
al. 1998a). Studies of HI lines in pre-main sequence stars 
have been extended also in the near and mid IR (Folha & 
Emerson 2001, Benedettini et al. 1998, Nisini et al. 1995, 
Natta et al. 1988). While the IR HI line ratios are well re- 
produced by wind models, they apparently fail to explain 
the observed profiles of the Br7 and Pa/3 lines. On the 
other hand, magnetospheric accretion models were able 
to reproduce the Br7 profiles showing redshifted absorp- 
tion components (MuzeroUe et al. 1998b) but not all the 
variety of profiles exhibited by T Tauri stars (Folha & 
Emerson 2001). Finally, very recently, spectro-astrometric 
observations of the Pa/3 line in four T Tauri stars, have 
indicated that the emission originates in the outflowing 
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material (Whelan et al. 2004). 

While there is ample literature on the origin and obser- 
vations of HI lines in T Tauri and Herbig stars, the inter- 
pretation of HI lines in more embedded Class I sources, 
for which optical observations of the Balmer series are 
impeded, is much less defined. Observations of spectrally- 
resolved HI lines in low-mass Class I sources so far have 
been limited to the Br7 line (Najita et al. 1996, Davis et al. 
2001). At variance with Ha in T Tauri stars, the profiles 
of the Br7 line appear more symmetric, with no evidence 
of the typical P Cygni signature. On this basis Najita et 
al. (1996) argue against a wind origin, while in the frame- 
work of magnetospheric accretion models, MuzeroUe et 
al. (1998b) use the Br7 luminosity in a sample of Class I 
sources to infer the mass accretion rates. 
To better constrain the emission region and the mech- 
anism for the excitation of HI lines in Class I sources 
it is necessary to obtain a larger observational set, in- 
cluding both information on line profiles and flux ratios 
among different lines. The aim of this work is to start 
such an analysis using both the ratios from lines in the 
Brackett series (the Brackett decrement) and high resolu- 
tion spectroscopy of different Brackett lines. The source 
investigated here is HHIOO IR, a low mass (Lhoi ^ 10 Lq, 
Wilking et al. 1992) Class I object located in the R CrA 
star forming core (D=130 pc, Marraco & Rydgren 1981). 
This source already has been recognized as an embedded 
protostar with high circumstellar activity on the basis of 
its near IR specrum (Greene & Lada 1996, Nisini et al. 
2004), and as also suggested by the fact that it is the 
driving source of an Herbig Haro object. It therefore rep- 
resents a sintable test case to perform a detailed analysis 
of its HI emission. 



2. Observations 

The observations were performed with the ISAAC spec- 
trometer, at the VLT UTl telescope, during 12 and 13 
July 2002. In the low resolution mode observations we 
employed a 0'.'6 slit to acquire spectra covering the H 
and K bands at R'^SOO. In the medium resolution mode, 
two spectral segments were acquired centered at 1.629/zm 
and at 2.161/Ljm and covering about 0.6/zm each. For 
these observations a O'.'S slit was used, giving a resolution 
of ~8000 and 9000 at 2.1 and 1.6/L(m, respectively. For 
both low and medium resolution observations, spectra 
of a standard star of B spectral type were obtained at 
an airmass similar to the scientific spectra, to correct 
for telluric absorption and obtain flux calibrations. The 
telluric spectra were carefully cleaned of any intrinsic 
HI absorption feature before being used. Wavelength 
calibrations were performed both using a Xenon lamp 
spectrum taken at the end of the night, and reflned each 
time on the OH sky lines observed in the spectra. This 
procedure leads to a wavelength calibration error of ~ 0.1 
A (i.e. about 1-2 kms~^). 



Table 1. Low resolution observations 
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"Fluxes and their errors are expressed in 10 erg cm ^ s 
Table 2. High resolution observations 
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" Vlsb. corrected for the velocity of the R CrA cloud (5.8 
kms~^ , Hariu et al. 1993) 



In the low resolution spectra we detected the Br7 at 
2.166/im in the K band spectrum, and lines from the 
higher levels of the Brackett scries, from Br 9 up to Br 
24, in the H band (Figure 1 and Table 1). The Br23 fine, 
which shows a line fiux much higher than the adjacent 
lines, is blended with a feature of Mg I. The medium res- 
olution spectra covered the Br7, the Br 12 (1.6113/Ltm) 
and the Br 13 (1.6411/xm) lines, which all appear resolved. 
Figure 2 shows the profiles of the lines, while in Table 2 we 
report the velocity information {Vlsr and AV (FWHM) 
) derived through a Gaussian fit. We estimate an error 
in the line width dc;termination that is between 5 and 
10km s~^ , based on the difference between the FWHM of 
the Gaussian fit with that directly measured with respect 
to the observed peak. The Vlsr velocities are corrected 
for the systemic velocity of the R CrA cloud taken as 5.8 
kms-i (Harju et al. 1993). 
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Fig. 1. Low resolution spectra in the H band of HHIOO IR, showing the emission lines from the Brackett series. 
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Fig. 2. Resolved spectra of the Br7, Br 12 and Br 13 lines 

in HHIOO IR. The imdcrlying continuum of each spectrum 
has been subtracted and each line is normalized to its peak 
intensity. The Vlsr is relative to the stellar velocity (i.e. 
the velocity of the associated molecular cloud has been 
subtracted). 



3. Line ratios 

The observations of a large number of HI lines from the 
Brackett series allow a detailed study of the decrement of 



this series, from which information on the line excitation 
can be derived. In normal HII regions, the decrement is 
actually used to derive the extinction through the emit- 
ting region, assuming that Case B recombination holds 
(Hummer & Storey 1987). However, the excitation of NIR 
HI lines in Young Stellar Object (YSO) envelopes usually 
cannot be reconciled with a standard Case B HII region 
since, due to the high densities of the source circumstel- 
lar gas, the lines may remain optically thick up to high-n 
levels. Therefore, to proceed with a correct analysis of the 
observed line ratios we first need to correct the observed 
lines for the reddening. Estimates on the visual extinc- 
tion for HHIOO IR, based on the optical depth of both 
the silicate feature at Q.l/xm and the 3/Ltm ice band ab- 
sorption feature, suggest a value around 25 mag (Whittet 
ct al. 1996). Such a value is not affected by the IR vari- 
ability of HHIOO IR (Molinari et al. 1994) since Graham 
(1998) found no variability in the 3/im ice band over a 
period of 7 yr. In addition, an between 20 and 30 mag 
has been also estimated from the comparison among the 
intrinsic and observed (H-K) colors, assuming the stel- 
lar spectral type derived from the photospheric features 
detected in the IR spectrum of the source (Nisini et al. 
2004). Therefore a value of Ay =25 mag, together with 
the reddening law by Rieke & Lebofsky (1985), have been 
adopted to correct the observed fluxes for the reddening. 
Figure 3 shows the ratios of the different Brackett lines 
with respect to the Br7 line intensity. In the same fig- 
ure the ratios expected for a Case B recombination are 
also plotted for the extreme cases of T=6000 K, n=10^ 
cm~"^ and T=10000 K, n=10^ cm"'', taken from Storey 
& Hummer (1995). The observed ratios are always higher 
than the Case B values, which indicates that the Brackett 
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lines remain optically thick up to high values of the n 
number. The displacement from the Case B curves could 
be due to a wrong extinction vahic assumed, but an Ay 
value as low as ~5 mag would be needed to reconcile the 
observations with the Case B recombination. 
In the figure, we also show the Brackett decrement in the 
case of emission from optically thick (r >>1) lines orig- 
inating from a region of fixed radius at T=10000 K. In 
this extreme case, lines from high n-number remain al- 
ways brighter than the Br7. It is evident from this figure 
that the observed ratios lie in an intermediate situation 
between these two extreme cases. Since for lines of a given 
series r oc A, the optical depth also decreases as the n- 
number increases. As a consequence, assuming that each 
line originates from a surface at which t '--^ 1 , it is expected 
that different lines trace a different emitting region whose 
size decreases as the n-number increases (see e.g. discus- 
sion in Benedettini et al. 1999). 

4. Line profiles 

The profile of the Br7 line in HHIOO IR is broad (FWHM 
^ 200 kms~^ ) and nearly symmetric. This profile is in- 
deed similar to those observed from other Class I sources 
(Najita et al. 1996, Davis et al. 2001) and from many T 
Tauri objects (Folha & Emerson 2001). These previous 
studies also showed that the Br7 line peaks are often 
blueshifted, as observed also in our case, where a Vlsr 
of —14 kms~^ has been measured. Such blue-shifted 
symmetrical profiles are difficult to reproduce either by 
wind models and by magnctosphcric accretion models. 
Wind models generally predict red-shifted peaks and 
P-Cygni absorption features. Most of the profile calcula- 
tions are however done for the Balmcr lines (e.g. Calvct 
et al. 1992), while predictions for NIR lines are seldom 
found. The Br7 and Bra profiles derived from Hartmann 
et al. (1990) in their magnetically driven wind model, do 
indeed lack P-Cygni absorption features and arc more 
symmetric, with a FWHM of the order of 200 kms^^ 
, comparable to the values observed in Class I sources, 
but they still present pronounced asymmetries and they 
are mostly red-shifted. In the case of wind models, a 
blue-shifted centroid may result from the occultation of 
the redshifted part of the profile by an optically thick 
circumstellar disk, as is observed in optical forbidden 
lines (e.g. Edwards et al. 1987). Such an effect however 
should be accompanied by a significant asymmetry in the 
redshifted side of the profile which apparently is lacking 
here. Magnetospheric accretion models, on the other 
hand, tend to predict profiles with blue-shifted peaks, as 
in our case, but where such a blue-shift centroid is the 
result of strong asymmetries caused by the suppression 
of the redshifted part of the profile (e.g. Muzerolle et al. 
1998a). Recently, Muzerolle et al. (2001) provided more 
refined line profile models, including line damping due 
to different broading mechanisms, which have the net 
effect of filling in the red-shifted absorption component 
causing a much more symmetric and centrally peaked 



profile. This effect is however particularly important 
only for the Ha line and becomes less significant for the 
higher Balmer and near IR lines. More symmetric profiles 
of the Br7 may be predicted at privileged inclination 
angles of the accretion disk (Muzerolle et al. 1998a), but 
even in the extreme situation of i '^70° the predicted 
red-shifted asymmetries should be still recognizable in 
our high S/N spectrum. Finally, the absorption of the 
red-shifted side of the profile also cause the FWHM line 
width to narrow significantly, and indeed the Br7 profiles 
shown in Muzerolle et al. (1998a, 2001) always have a 
FWHM of about 100 kms"-'^ or less, thus much narrower 
than the values we observe. Such a discrepancy between 
the FWHM predicted by the magnetospheric accretion 
models and the observed wider Br7 line-widths was also 
pointed out by Folha & Emerson (2001) for their sample 
of T Tauri observations. Since the line width depends on 
the infall velocity field in the accretion flow, which in 
turn is a function of the source gravitational potential 
energy, such a width could be larger in sources with a 
mass greater than the canonical value assumed by the 
Muzerolle et al. models, i.e. 0.5 Mq. In the specific 
case of HHIOO IR, a stellar mass of 0.4 Mq has been 
estimated (Nisini et al. 2004), so in principle line widths 
even narrower than those predicted by Muzerolle et al. 
should be expected. 

On the basis of these considerations, we believe that at 
present both the existing wind and accretion models still 
fail in reproducing the observed HI near-IR symmetric 
profiles, and therefore it is not possible to favor one of 
these two excitation mechanisms on the basis of profiles 
alone. 

The profiles of the observed Br 12 and 13 lines also ap- 
pear broad but with FWHM velocities smaller than the 
Br7 linc\ These profiles are also rather symmetric: two 
apparently blue-shifted absorptions in the Br 13 line can 
be ascribed to the presence of two photospheric absorption 
features, as indicated in Figure 3. If we correct for these 
absorption features, the Brl3 line-width becomes ~190 
kms^^ , i.e. more similar to the Brl2 line and still sig- 
nificantly lower than the Br7 line-width. It seems there- 
fore to be a trend for the FWHM velocity to decrease 
with higher quantum numbers. We have seen that the 
observed Brackett lines are thick and have different op- 
tical depths, thus they trace zones at different physical 
depths in the emitting region. Since moreover the opti- 
cal depth decreases with increasing upper quantum num- 
ber (i.e. with decreasing wavelength), it is expected that 
high-n lines trace regions more internal than the Br7 line. 
Therefore if we assume that the observed line widths mea- 
sure the maximum velocity attained at the emitting region 
surface, the observed trend would imply that the velocity 
is increasing going outwards, as expected, e.g., in an ac- 
celerating wind. With the assumption that such a veloc- 
ity trend can be ascribed to the gas kinematical motion 
alone, then the different observed linewidths are incom- 
patible with a velocity law due to any accretion process or 
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Keplerian rotation (such as ii oc r~^/^) where the velocity 
is decreasing with the distance from the central object. 
Emission from a flattened geometry viewed at different 
inclination angles cannot produce changes in this trend. 
Muzerolle et al. (2001), on the other hand, presented ev- 
idence that other line broadening effects, radiative and 
Stark broadening in particular, may dominate the line- 
width in some circumstances, an effect which is particu- 
larly important for the Ha line. However. The Muzerolle 
et al. (2001) models show that these broadening effects 
become negligible for near IR lines, and hence the line- 
widths of the Brackett lines arc more likely to reflect the 
dynamical properties of the gas in the emitting region than 
other broadening effects. 

A rough estimate of the size of the emitting region can 
be given assuming that the flux in the different lines is 
given by a blackbody emission at T~10000 K, i.e.: 



JJ-^ c 



(1) 



Taking as AV' the widths of the observed lines and a dis- 
tance of 130 pc, we derive that the Br7 emission size is 
of the order of 4 10^^ cm while the Brl2,13 emission sizes 
are of the order of 2 10^^ cm. This indicates that the HI 
lines originate from a region of only 4-6 R0, thus very 
close to the protostellar photosphere, which is expected 
to be a few solar radii in diameter. Since the emitting re- 
gion is not vcy large, the central source may produce an 
occultation of the redshifted part of the flow. A simple 
way to check the amount of this occultation is through 
geometrical considerations. The emission originates from 
the surface at which r '--^1, thus the amount of flux not 
reaching the observer due to the central star can be es- 
timated as the ratio between the emission region surface 
and the source projected area. This estimate shows that 
only ^ 5 % of the total line emission would be hidden by 
the central source, and thus difficult to detect although it 
should mainly affect the line redshifted wing. An occul- 
tation fr(mi a circumstellar disk should in fact produce a 
more severe effect on the line profile. Such an effect is not 
observed also in HI lines from T Tauri stars while it is 
often present in other optical forbidden lines such as [01] 
(e.g. Edwards et al. 1987, Takami et al. 2001). This has 
been interpreted as evidence for the presence of internal 
disk holes allowing the redshifted part of lines coming from 
compact regions to reach the observer while obscuring that 
of the forbidden lines, originating from more external and 
lower density regions (Takami et al. 2001). 

5. Comparison with models of ionized envelopes 

There are very few models in the literature treating in 
detail (i.e. with a full radiative transfer treatment and 
taking into account geometrical effects) the excitation of 
HI IR lines and providing predictions for both flux ratios 
and line profiles. Natta et al. (1998) provide Pa/3, Bra 
and Br7 intensities of T Tauri stars in the framework of a 
partially ionized wind model but without any prediction 



of line profiles. Muzerolle et al. (1998a), on the other 
hand, compute the Br7 profile expected from their 
magnetospheric accretion model, but the line fluxes are 
derived only for Pa/3 and Br7 lines. 

Giving the difficulty of comparing our data with specific 
models, wc have adopted a more general approach, 
considering a simple model of an ionized envelope in 
which the gas velocity follows a general law of the type 
V = Vq + {Vmax — ^)(1 — (^i J which assumes 
that the gas is accelerated at a maximum velocity Vmax 
at a distance that depends on the parameter a. The 
ionized envelope extends from an initial radius to an 
outer radius tq. Inside this region, the electron density 
follows from the continuity equation and it is equal to 
rie = Mxe/47rmjy y(r), where is the fractional 
ionization and M is the rate of mass flow inside the 
region. The radiation transfer has been treated in the 
Sobolev approximation assuming the gas in LTE and 
following the formalism described in Nisini et al. (1995), 
where a discussion about the limitations of the adopted 
assumptions is also given. 

We have tuned the input parameters in such a way as 

to reproduce the line ratios, line fluxes and the observed 
line widths with the assumption that the FWHM traces 
the maximum velocity attained in the line emitting re- 
gion. This last constraint is equivalent to reproducing the 
emitting region sizes as estimated in the previous section. 

In this framework, the optical depth at line center, 
which is the main parameter affecting the line decrement, 
is given, as a function of the distance from the initial ra- 
dius n, by: 



MXe 



V{r)- 



(2) 



where is the line absorption coefficient and p the line 
center frequency. From this expression, we see that the 
optical depth is strongly affected, by the adopted velocity 
law, the initial radius and the rate of ionized gas flow 
Mion=M Xe- We have chosen that the emitting region be- 
gins at the protostellar radius _R* . assumed equal to 3 Rq . 
This is the predicted radius of protostellar objects and it 
is consistent with the derived value of (3.5±0.7)Rq esti- 
mated from the inferred spectral type and stellar luminos- 
ity (Nisini et al. 2004). 

With this condition, to maintain the lines optically 
thick up to high-n values, the Mjon value needs to be suf- 
ficently high, e.g. larger than 10~^ Mq yr"-*^. At the same 
time, the absolute line fluxes are also sensitive, in addition 
to Mion; to the gas temperature and the emitting region 
size. Assuming T= 8000 K, we find that ro values of a few 
stellar radii are needed to mantain the observed line fluxes 
close to the observed absolute values. Colder emission may 
be consistent with slightly larger emission regions, which 
however cannot exceed ~8 i?* for T ^ 4000 K. 

In Fig. 3 we show the predicted decrement in our best 
fit model with the following parameters: M Xe=2 10~^ 
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Mgyr-i, ro=3R*, Vo=30 kins"! , V„a^=230 kms"! 
and a=4. Such a model also reproduces fairly well 
the dereddencd flux of the Br7 line (e.g. 5.9 10~^^ 
ergcm~^s~^ vs (6.2±.2) 10""'^^ ergcm"^ s"""^); moreover it 
predicts that the maximum velocities corresponding to the 
radii where the Br7, Brl2 and Brl3 lines attain t=1 arc 
220, 190 and 180 kms^^ respectively, also in agreement 
with the observed line widths. The model however overes- 
timates the emission from high Hup numbers with respect 
to the observed values. The high n^p lines originate from 
inner regions where, according to the adopted velocity law, 
the speed is still low; consequently, the expected optical 
depth of these lines is mantained high despite their higher 
frequency, according to the relationship (2). A better fit to 
the data can be obtained by assuming a different velocity 
law with a much higher initial velocity (~ 180 kms~^ ), 
which however seems unphysical implying that the ionized 
gas is flowing from the stellar surface already at a velocity 
close to the escape speed. 

In summary, and considering the oversemplification of 
the adopted model, we can reasonably conclude that the 
different observing features can be accounted for by an 
ionized gas, emitted in a region of a few stellar radii close 
to the source where the gas is already accelerated at about 
200 kms~^ and it is expanding, after an initial steep ac- 
celeration, with a rate of (ionized) mass loss of the order 
of 10-"^ MQyr-i. 

In deriving this mass loss rate value, it is assumed that the 
ionization fraction does not change over the line emitting 
region. While this may not be adequate if one considers 
the wind ionization structure due to the stellar irradiation 
alone (e.g. Natta et al. 1988), it is a rather acceptable ap- 
proximation in the present case, given the fact that the 
lines originate in a compact emitting region close to the 
source, and that other heating mechanisms have been rec- 
ognized to be important in maintaining a high degree of 
ionization in coUimated winds (e.g. Shang et al. 2002). 

The origin of such an emission in a very compact region 
close to the source suggests excitation from a stellar wind, 
more than from a disk-wind unless the emission comes 
from the inner section of a disk extending down to the 
source surface. This evidence naturally excludes an origin 
in MHD disk-winds (e.g. Ferreira & Pelletier 1993) where 
the disk is truncated at large distances from the source by 
the action of a strong magnetic field. 

The Mion value derived in our model is moreover large 
enough to suggest that the gas ionization in the emit- 
ting region should be high. An ionization fraction of ~ 0.1 
would imply an M ~ IO^^Mq yr~^. Values of the order 
of 10^^-10~* Mq yr"-'^ have been estimated in T Tauri 
stars (see e.g. Natta et al. 1988), it is however expected 
that the efficiency of the mass loss mechanism is higher 
for younger sources. 

Finally the adopted velocity law impies that in the out- 
flowing gas the velocity rapidly increases with the distance 
from the central source and it is already high (~200kms^^ 
) at about two stellar radii. This velocity trend, which 
is consistent with magnetically-driven stellar winds (e.g. 



Lago 1984, Hartmann 1990) as well as with magnetocen- 
trifugal wind models like X-winds (Shu et al. 1994) has 

been also measured in some T Tauri stars from spectro- 
astrometric observations of the emission (Takami 2001, 
2003). 

Our simplified model, in which LTE emission is 
assumed, necessarily predicts symmetric and centrally 
peaked emergent line profiles. The P Cygni profiles and 
blue-shifted absorption components which seem to be a 
common characteristics of wind models thus far explored, 
result when large line opacities are coupled with significant 
deviations from equilibrium in the external expanding re- 
gion of the wind, where the increasing inefficiency of colli- 
sions causes a decrease of the source function. On the other 
hand, the mere fact that the observed lines appear sym- 
metric and almost Gaussian could be an indication that 
LTE conditions for the Brackett lines actually hold in the 
region where they originate. Indeed, Nisini et al. (1995) 
compared their LTE model with models treating in detail 
the hydrogen level populations, reaching the conclusions 
that the LTE assumption is a reasonably good approxi- 
mation of the Brackett lines, while it may be not totally 
adequate for the Pa/3 line. Moreover, the limited Br7 line 
profile predictions presented by Hartmann (1990) shows 
that since the upper hydrogen levels reach equilibrium 
more easily than the Balmer lines, the resulting absorp- 
tion features are much weaker and the degree of asymme- 
try significantly reduced. Obviously a more refined model 
than the one presented here, with a full statistical equilib- 
rium treatment and a more realistic geometry, is needed 
to better explore the ability of wind models to reproduce 
the observed Brackett symmetric profiles. 

6. Conclusion 

The analysis of the Brackett decrement in the HHIOO IR 
Class I source has shown that the observed Brackett lines 
remain optically thick up to a high quantum number. In 
addition, the profiles of the Br7, Br 12 and Br 13 lines are 
fairly symmetric with a tendency of the line FWHM to 
decrease as the n-mimbcr increases. On the basis of this 
evidence we argue that a wind-like mechanism, where the 
gas is accelerating outwards is more suited to reproduce 
the observed features. A comparison with a very simple 
model suggests that the emission region should be very 
compact and close to the stellar surface, a fact favour- 
ing a stellar wind or the inner region of a disk wind with 
a small truncation radius as the natural emission region. 
The mass loss rate of the ionized gas should be fairly high, 
of the order of Mion ~ 10~^Mq yr"-'^, to maintain the lines 
optically thick. This suggests that the ionization fraction 
in the emitting gas should not be smaller that ~0.1 to 
have a total mass loss rate not exceeding typical values 
expected from Class I objects. 

Our analysis, however, does not solve the problem of how 
to reproduce the observed symmetrical line profiles of op- 
tically thick IR lines, which at present are apparently dif- 
ficult to modell by both wind and accretion models. This 
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probably points to the fact that the real situation is more 
complicated than described in the simple model presented 
here. Another caveat to the derived conclusions is the near 
IR photometric variability of HHIOO IR, which shows dif- 
ferences in the K band of more than a magnitude over a 
period of a few years. Such a variability can also affect the 
line emission of this source. Indeed, variability not only 
of HI NIR line fluxes but also of their profiles has been 
observed in other Class I sources exhibiting photometric 
variability (Nisini et al. 1994). 

Finally, our conclusions are valid for the specific case 
of HHIOOIR and as such cannot be generalized to other 
Class I objects. Indeed, more asymmetric profiles of the 
Br7 line in similar sources have been succesfuUy modelled 
by magnetospheric accretion models, such as e.g. WL16 
(MuzeroUe et al. 1998b). For a better understanding of 
the origin of the HI lines in Class I sources, a more de- 
tailed modelling effort is needed, allowing us to fit in a 
consistent way both the intensity and profiles of the dif- 
ferent observed lines. 
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